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ABSTRACT 


A three-dimensional version of the direct simulation Monte Carlo method is adapted to assess the con- 
tamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite. 
Emphasis is placed on simulating a realistic, worst-case set of flowfield and surface conditions and geometric 
orientations for the satellite in order to estimate an upper limit for the cumulative level of volatile organic 
molecular deposits at the aperture of the Halogen Occultation Experiment (HALOE). Results pertaining to 
satellite environment are presented regarding contaminant cloud structure, cloud composition, and statistics 
of simulated molecules impinging on the HALOE aperture, along with data related to code performance. 
Using procedures developed in standard contamination analyses, along with many worst-case assumptions, 
the cumulative upper-limit level of volatile organic deposits on the HALOE aperture over the instrument’s 
35-month nominal data collection period is estimated at about 13 , 350 A. 
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1.0 INTRODUCTION 

The Upper Atmosphere Research Satellite (UARS) was designed to collect data on various aspects of 
Earth’s upper atmosphere to help characterize global atmospheric changes that are thought to be occurring. 
Figure 1 depicts the UARS model used in this investigation. One instrument carried on board the currently- 
operational satellite is the NASA Langley Halogen Occultation Experiment (HALOE). During each sunrise 
and sunset experienced by the satellite, the HALOE instrument peers at the Sun through the atmosphere 
and measures solar infrared energy absorption over a range of wavelengths. With this information, vertical 
distributions for concentrations of trace gases such as HF , HCl , and 0 3 are deduced. 1 

In order to obtain precise measurements, HALOE optical surfaces must be kept extremely clean. Experi- 
ence shows that even minute amounts of contamination can seriously degrade the spectral transmissivity and 
reflectivity of optical surfaces. 2 Although atmospheric density is extremely low at the satellite orbit altitude 
of 600£m, there are many potential sources of contamination contained on board UARS itself, ranging from 
light gases to heavy, long-chain polymers. Volatiles outgassed from certain materials used in satellite con- 
struction, gases vented from the interior of the satellite and cryogenically-cooled experiments, and effluents 
exhausted from attitude thrusters, all contribute to a gas cloud that envelops the vehicle. 

Previous analyses have indicated that the cumulative level of deposits on critical HALOE optical sur- 
faces over its nominal 35 months of operation could exceed maximum acceptable levels by over an order 
of magnitude. 3 While UARS instrument managers specified a goal of no more than 100 A of cumulative 
deposits at the HALOE aperture, analyses predicted for its original set of orientations that it would receive 
about 32,170A. 3 A recommended change adopted in stowing positions for HALOE reduced that estimate to 
about 6 188 A , 3 These analyses utilized SPACE2 4 and MOLFLUX, 5 which are standard engineering codes 
developed for estimating fluxes on target surfaces due to outgassing. These codes either assume no inter- 
molecular collisions within the gas cloud or employ first-collision theory. 

Because of the important effect of contamination on the success of the HALOE data collection effort, 
detailed characteristics of the gas cloud in the vicinity of the HALOE aperture are needed to accurately 
estimate incident contaminant fluxes and to help define criteria for optimal operation of the instrument. In 
order to facilitate this goal, a study of the contamination environment surrounding UARS was performed 
using the direct simulation Monte Carlo (DSMC) technique. 6,7 

This report describes details of the investigation, emphasizing methods used and assessment of input data 
and assumptions. First, a discussion of issues important to satellite contamination is presented, followed by a 
general description of the DSMC method and preprocessors specifically developed to handle various necessary 
input conditions and geometric orientations required for the investigation. This discussion is then followed 
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by a description of specific data pertinent to contamination of the HALOE instrument on board UARS. 
Results include simulated flowfield species number density maps describing the overall gaseous environment 
surrounding UARS, and statistical details regarding the fluxes intercepted at the HALOE aperture. Also, 
results regarding predictions of cumulative contaminant deposit levels are presented and interpreted, and 
code performance statistics are mentioned. 

2.0 CONTAMINATION ISSUES 

Satellites tend to surround themselves with their own artificial atmosphere through surface outgassing, 
equipment venting, and attitude thruster firing. This artificial atmosphere typically has very low density, and 
its overall distribution is usually anisotropic and non-Max wellian. It is composed of species having a wide 
range of molecular weights, including massive, long-chain, volatile organic compounds, and is characterized 
by a wide range of concentrations, with some important species present only at trace levels. 

The molecular flux of contaminants incident on a spacecraft surface has traditionally been subdivided into 
two categories: direct and return flux. Freestream or outgassing molecules directly impinging on a critical 
target surface that may have encountered prior collisions with other surfaces, but no intermolecular collisions 
with other gas species, belong to the direct flux. Outgassing molecules having undergone intermolecular 
collisions before reaching the target surface belong to the return flux. Direct flux computation is relatively 
straightforward for diffuse surfaces (geometric vie wf actors or particle tracing), but return flux computation, 
which requires an integrated knowledge about the contaminant cloud structure and collision physics, is 
intrinsically more complex. It also requires that a relatively large computational domain surround the 
satellite to ensure that the simulation satisfactorily accounts for virtually all molecules which can return to 
the vehicle surface, regardless of how far they ventured away from the satellite before being scattered back. 

The physics involved in contamination processes, such as surface desorption, reflection, absorption, ad- 
sorption, etc. are very complex. The experimental database is presently insufficient to devise accurate models 
to simulate mechanisms such as surface outgassing (flux magnitude as a function of surface temperature and 
material composition, emitted angular velocity distribution, and species composition), surface deposition 
and polymerization (effects of ultraviolet radiation, atmospheric monatomic oxygen flux, and surface mate- 
rial), and surface reflection (momentum and energy accommodation, reflected angular velocity distribution). 
Reliable estimates of the contamination environment surrounding future spacecraft will therefore require 
improving the present state of knowledge. The molecule- based DSMC algorithm is readily amenable to 
incorporate new physical models as these become available. DSMC could therefore be used to evaluate these 
models and foster new research in this field. 
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3.0 NUMERICAL TOOLS 

In this section, the method used to investigate HALOE on-orbit instrument contamination is presented, 
along with preprocessing routines developed to assimilate data concerning geometric configurations and 
outgassing required to provide detailed simulations. First, a general description of the DSMC method is 
presented, followed by a brief presentation describing the various preprocessing routines used, many of which 
have been specifically developed to handle outgassing emission data. 

3.1 DSMC METHOD 


The direct simulation Monte Carlo method has been steadily gaining acceptance as an approach to 
solving problems related to fluid dynamics when the effects of rarefaction become important and the gaseous 
medium no longer behaves as a continuum . 8 In contrast to the approach of Computational Fluid Dynamics 
(CFD), where the gas is treated as a continuum mathematical fluid, macroscopic flowfield properties are 
calculated based on the interactions of many thousands or millions of particles that behave as molecules at 
the microscopic level . 8 

In the DSMC method, the simulated flowfield volume is discretized into a network of subcells, and the 
ensemble of particles constituting the gas is allowed to develop in time, where time is discretized into small 
increments over which the processes of molecular collisions and molecular motion are decoupled. The subcells 
are used to help identify nearest-neighbor prospective collision partners. The flowfield evolves in time from an 
initial, specified configuration, and once it reaches steady state, cumulative statistics on the gas are kept by 
the computer for small numbers of subcells grouped together to form cells. By specifying the intermolecular 
potential governing the behavior of collisions between the simulated molecules, or “particles,” and by using 
concepts of kinetic theory and statistical mechanics to sample the flowfield at the subcell or cell level, one 
can establish macroscopic flowfield characteristics. 

Although the true intermolecular potential for most interactions are unknown, a number of simple models 
are available which give reasonably accurate macroscopic behavior. The potential commonly used in con- 
temporary DSMC applications is the Variable Hard Sphere (VHS) model . 8 In this model, the hard-sphere 
assumption of uniform, isotropic scattering is retained, but collision cross-sections become proportional to 
the relative velocities between molecules . 8 The VHS model has been employed in this investigation. 

The DSMC outgassing model has been widely used in other applications , 9 ' 10 and is assumed to be valid in 
this investigation. It assumes that molecules of an outgassing species are emitted from specified surfaces at 
given rates, with velocities sampled from a drifting Maxwellian velocity distribution at a given temperature. 

The DSMC code used in the present study was devised by Bird 11 and modified by Rault . 6,7,12 High com- 
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putational efficiency is achieved with this code through the use of an unstructured cell network overlaid on 
a cubic Cartesian sub cell mesh. Two levels of subcells are employed in this particular code. A fine scale 
mesh (FCM) is used to resolve high fiowfield gradients that often occur near windward-side surfaces, and a 
coarser mesh (CCM) is used farther away to conserve computational resources and to improve statistically- 
based calculations in regions where overall density levels are low. The individual subcells at the FCM level 
are referred to as “pixels.” When available, comparisons between results generated with this code and exper- 
imental measurements regarding surface and fiowfield properties for flows around a wide variety of vehicles 
have yielded excellent agreement. 12 - 13 The code has recently been optimized for vector architecture super- 
computers, and has been complemented with a set of utilities for diagnosis, preprocessing, postprocessing, 
and grid adaptation. 12 Grid adaptation of the cell network with respect to certain variables can enhance 
resolution of fiowfield phenomena as well as surface fluxes. 

Special coding has been developed to collect information regarding particles striking the H ALOE aperture. 
Such data include particle velocity, position where the aperture plane has been struck, species type, etc. In 
order to distinguish molecules belonging to the direct flux or the return flux, a parameter has been added 
to the run code to store the number of surface and intermolecular collisions each particle encounters prior 
to reaching the aperture. An additional parameter, storing information identifying the originating surfaces 
of these outgassed molecules, helps indicate which surfaces dominated the influence of volatile species fluxes 
at the aperture plane. 

A novel additional capability featured in this simulation is the introduction of a parallel processing tech- 
nique to separate the fiowfield into two domains. This option allows higher resolution in the inner domain, 
combined with a larger outer domain for computing return flux, than would be possible if the resources of 
only a single computer were available. Implementation of this technique is discussed further in Section 4.2. 

3.2 PREPROCESSOR ROUTINES 

3.2.1 Geometry Definition 

To define and input the inherently complex geometry of typical satellites, the interactive three-dimensional 
CAD preprocessor developed by Rault is used.®’ 7 This preprocessing utility relies on the fact that even a very 
complex body geometry can be decomposed into a series of simple geometric primitives, such as spheres, 
cylinders, planes, etc. or portions thereof. 6 - 7 The current set of 16 basic primitives is presented in Table 1. 
UARS is modeled using 385 of these primitive subelements, and the resulting geometric model is displayed 
in Fig. 1. Also shown in the figure is the relative coordinate system used in this investigation, and details of 
the HALOE instrument in Figs, lb- Id. As described in Ref. 6, properties of these subelements are entered 
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into the preprocessor code through an input file. Typical input data are presented in Table 2. After listing 
data associated with the overall configuration of the flowfield volume, a series of fields describe the type of 
each subelement, its dimensions, location, and orientation. Since the DSMC code assumes SI units, use of 
any other system requires conversion using the “scale” factor. In Table 2, lengths had been measured and 
entered in inches, requiring “scale” to be set to 0.0254. 

Other entries displayed in Table 2 also require further explanation. Certain types of geometric elements 
can be completely described by quadratic functions. These elements are categorized under “quadrics,” and 
are handled analytically by the DSMC code. Other planar surfaces, such a s rings and disks, where limits 
cannot adequately be described by quadrics, are simply referred to as “planes,” and are stored as a number 
of points with high resolution. This resolution must be higher than that used to discretize the surfaces 
for incorporation into the DSMC code for high fidelity. 1112 In this study, these two groups encompassed all 
types of subelements used. In other studies, nonplanar surfaces of arbitrary complexity are defined in similar 
fashion to the rings and disks used here. 12 

Since subassemblies or components constituting the satellite are represented in the code by groups of 
geometric subelements, these subelements may be linked by a common “component number.” Defining a 
unique color to represent each component number can be useful in geometry visualization to highlight or 
otherwise distinguish individual components. References to “limit planes” in Table 2 pertain to the fact that 
these surfaces are discretized within the DSMC framework and replaced by small cubes. 11 Along certain 
edges where surfaces should intersect, small gaps may occur as a result of the discretization process. This 
problem occurs most frequently where surfaces intersect at oblique angles. Since the flowfield cell network 
is built up in layers from surfaces, these gaps will allow the interior to be filled with flowfield cells, which 
consumes computational resources and may lead to confusing errors. 12 Limit plane displacement lengthens 
or widens the original subelement to eliminate those gaps. Finally, the entry fields identified as “F3I limits” 
refer to a subordinate volume contained within the flowfield where the particular subelement lies. This 
information is used to speed up the discretization process by telling the code to check for a particular surface 
only within that subvolume instead of searching the entire flowfield. 

Dimensions and locations for various facets of the UARS satellite have been obtained from engineering 
drawings provided by HALOE instrument managers, and a complete listing of the data entered for the 
baseline geometric representation is presented in Table 3. Code has been developed to alter subassemblies of 
the geometric UARS model by taking the baseline geometry file listed in Table 3, and rotating the required 
subelements about a given axis at a given point. The input files for cases studied are displayed in Table 4, 
and the output resembles Table 3, except that values for the affected subelements have been duly altered. 
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3.2.2 Out gassing Inputs 

Code has been developed to process information pertaining to individual outgassing surfaces for incorpo- 
ration into the general DSMC run code. These input files include subelement-level information pertaining 
to the outgassing flux at a given reference temperature, and the molar composition of the outgassing flux. 
Table 5 is an example of typical outgassing information, and Table 6 is a reproduction of the input file for 
outgassing information used in this investigation. The values listed before each outgassing rate are species 
mole fractions for each of the six species modeled. 

The preprocessor code contains enough flexibility to process information on whether or not an individual 
subelement is composed of several regions outgassing different species at different rates. This capability has 
been used to create vent sources on certain surfaces of the Instrument Module (IM) by superimposing vent 
locations, flux rates, and gas constituents on existing larger elements. An example is illustrated in the data 
fields associated with the first geometric element in Table 6. Since outgassing rates are strongly dependent on 
surface temperature, additional coding has been developed to assign a unique value to each satellite surface. 
A listing of subelement-level temperatures in °C is presented in Table 7. 

3.2.3 Vent Plume Modeling 

In addition to outgassing information, a preprocessing routine for handling jet-like vent plumes from sonic 
orifices has been developed. Whereas effluents outgassed from surfaces and exiting most vents diffuse into the 
flowfield with modeled velocity distributions resembling drifting Maxwellian distributions at 7; urface , “jet” 
sources represent relatively high mass flux rates exiting small orifices as choked, uniform flow at a supplied 
exit temperature T ex it* The data necessary to recreate such sources are shown in Table 8. 

3.2.4 Global Resolution 


Finally, Table 9 lists inputs used by the DSMC code in order to initialize the overall simulation. Some 
of the input variables are simply switches used to tailor the code to the problem at hand with respect to 
where high grid resolution is needed, sizing maximum lengths for Fortran arrays, information required for 
setting reasonable estimates for DSMC simulation parameters, flow direction with respect to the body, and 
the like. 12 
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4.0 APPLICATION TO UARS/HALOE CONTAMINATION 

4.1 UARS/HALQE DESCRIPTION 

UARS was placed in orbit in September 1991 to probe Earth’s upper atmosphere and characterize global 
atmospheric changes which are thought to be taking place. It operates at a roughly circular 600ibm altitude 
orbit with a 96-minute period (r). Depicted in Fig. 1, UARS is about 10.3m long, 5.9m high, 3.2m wide, 
and its solar panel (SP) measures 13.3m by 3.3m. 

Several instruments have been mounted on the Instrument Module (IM), including HALOE, (an occul- 
tation viewing radiometer that requires part per million signal precision and characterization 1 ), the Cryo- 
genic Limb Array Etalon Spectrometer (CLAES), and the Nadir and Zenith Energetic Particle Spectrome- 
ters (NEPS and ZEPS). The High Gain Antenna (HGA) receives signals from UARS mission controllers, and 
the Solar/Stellar Positioning Platform (SSPP) assists in navigating the satellite. The Multi-mission Mod- 
ular Spacecraft (MMS) is used to periodically boost or maneuver the satellite in accordance with mission 
guidelines. 

UARS is designed to always keep the side with HALOE and the SP facing the Sun when it is over the 
horizon, but HALOE only takes measurements during sunrise and sunset, which occur once every period. 
At other times, whether the Sun is above or below the horizon, HALOE is stowed facing away from the 
freestream direction and the bulk of the satellite. Even though HALOE’s stowed position has been chosen to 
minimize contamination from surface outgassing, it will periodically have substantial viewfactors from the 
continuously-rotating solar panel. 

4.2 PARALLEL CODE SCHEME 


As described in Section 3.1, the flowfield has been separated into two domains in order to increase grid 
resolution near the satellite and to create a large outer domain for computing return flux. Figure 2 shows 
a schematic representation of the computing arrangement. A different processor computes the solution in 
each domain independently, and Parallel Virtual Machine software from Oak Ridge National Laboratories 
has been implemented to allow the processors to communicate data between one another. 14 

The potential utility of this parallel routine software has been demonstrated even though an outer domain 
was not considered necessary for the scope of this particular investigation. Although the effect of return 
flux on HALOE is deemed insignificant, this feature could be useful in other contamination studies (at lower 
altitudes, for instance), where the far-field effect could play a more significant role. In this investigation, the 
inner domain takes on the shape of a right parallelepiped with dimensions slightly larger than the maximum 
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dimensions of UARS, and the outer domain is a cube. During a given timestep At, particles reaching or 
passing through the interface boundary separating the domains are moved to the boundary, and information 
required by the other processor for these entering particles is stored in a buffer. Data include particle position 
and velocity information, remaining fraction of At, individual species codes, and individual tag parameters. 
The buffer file is transmitted as a message to the other processor once each timestep. 

It is reasoned advantageous to size the domains using the following equation: 

FNUM, _ At, Vi /a 3 

FNUMo A t 0 ~ V 0 “ [lj ' W 

Concentrating on the first equality in Eq. (1), the need to clone or remove particles at the interfaces to 
ensure correct transmitted species flux levels is obviated. The second equality is an attempt to achieve 
roughly equal numbers of particles in each domain, a desirable feature which could result in roughly equal 
processor loadings under certain circumstances. One compromise inherent in Eq. (1) requires finding suitable 
timesteps for the different domains, such that when adequate temporal resolution in the inner domain is 
obtained, particles do not travel too quickly through the outer domain. Of course, shorter timesteps translate 
into longer run times. 

4.3 SPECIES WEIGHTING FACTORS 

Due to the large variation in species concentrations within a contaminant cloud, “species weighting fac- 
tors” (SWF) have been implemented to artificially enhance minor species concentrations. Important species, 
such as the organic volatiles most likely to deposit and remain on critical surfaces, are usually present in trace 
amounts, and statistical sample sizes concerning these molecules would therefore be inadequate without the 
use of these weighting factors. The DSMC code has been modified by artificially increasing the outgassing 
flux of trace species by a species weighting factor, which in effect enables one to specify species-dependent 
values of FNUM to obtain reasonable numbers of each representative simulated species within the flowfield 
for better statistical property measurement. The only aspects of the simulation affected by this artifice 
are collision rates and quantities associated with flowfield sampling. These aspects may be accounted for by 
reducing bimolecular collision probabilities and cumulative species sampling summations by the same factors 
within the code. 

The conditions under which SWF’s may be applied should be considered carefully, since mean free path 
lengths of trace species may be artificially increased when intermolecular collisions become important. Such 
an approximation is not considered a problem for the simulated UARS/HALOE environment. 
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5,0 COMPUTATIONAL CONDITIONS 

This section describes details concerning the application of the tools mentioned above to the analysis of 
the UARS/HALOE contamination environment. First is a description of the geometric configurations that 
were considered to provide realistic worst-case estimates for deposits at the HALOE aperture, followed by 
assumed freestream and outgassing conditions. 

5.1 GEOMETRIC CONSIDERATIONS 

5.1.1 Configurations 

For the set of assumptions associated with this approach, a realistic, conservative estimate for the cumula- 
tive level of contaminants deposited at the HALOE aperture must recognize that certain instruments rotate 
with respect to the rest of UARS over time. The solar panel (SP) continuously tracks the Sun, and the 
Solar/Stellar Positioning Platform (SSPP) undergoes limited movement as it alternately tracks the Sun and 
certain stars. HALOE faces the Sun during sunrise and sunset every 96-minute orbital period. Otherwise, it 
is stowed facing away from the freestream direction with severely restricted views of the rest of the satellite. 3 
Three different configurations have been chosen, representing a realistic worst case set of orientations. Two of 
these represent the beginning and end of the HALOE data collection period considered most deleterious with 
respect to contamination, and one represents the instrument in its stowed position after this data collection 
period. 

Since the direct flux of outgassed molecules intercepted by a target surface is strongly dependent on 
the solid angle Q subtended by the originating surfaces within the target’s line of sight (viewfactor), cases 
have been chosen where HALOE collects data facing the -X-direction (toward the Multi-mission Modular 
Spacecraft (MMS)). In this orientation, viewfactors subtended by the SP and the bulk of the satellite are 
maximized. In addition, since both the SP and HALOE track the Sun during HALOE data collection 
periods, the SP presents a large frontal area to the oncoming freestream, which could potentially contribute 
to increased overall number flux levels at the HALOE aperture. Movement of the SSPP has been ignored. 

The angle between the HALOE aperture normal and the UARS body axis during data collection periods 
is equal to /?, the angle between the Earth-Sun line and the plane of the satellite orbit. 3 Angle /? also plays 
a role in outgassing, since satellite surface temperatures are directly proportional to local surface inclination 
with respect to the Sun. For this particular orbit (0° < f3 < 80°), temperatures of influential surfaces are 
generally maximized when /? = 80°. One important exception is that trends for the SP surface temperatures 
run counter to the other objects. Due to the glancing angle it makes with the Sun at 0 = 80°, SP surface 
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temperatures, and hence its outgassing rates, are minimized there. Conversely, SP surface temperatures are 
maximized when /? = 0°. The large SP surface area and its orientation with respect to HALOE make it a 
formidable influence on the latter’s contamination environment. Another influential factor is that although 
the temperatures of the many influential outgassing surfaces are minimized at = 0°, their viewfactors are 
maximized with respect to HALOE for this orientation. 

The three sets of geometric orientations studied in this investigation are listed below in expected descend- 
ing order of severity. Angle a represents azimuth angle with respect to the UARS X axis, and 0 represents 
elevation. Case 1 assumed a = -180° and 0 = 6° (end of sunrise data collection period); Case 2 assumed a 
= —180° and 0 = -23.6° (beginning of sunrise data collection period); and Case 3 assumed a ~ 45° and 0 
= -25° (stow position immediately prior to /? = 0° sunset data collection period). Views of the geometric 
configurations for these cases are depicted in Fig. 1. Note that the SP location in Fig. 1c lies behind the 
HALOE aperture plane. The visual absence of that large, relatively high-temperature, outgassing object 
nearby should have a significant impact on the Case 2 contamination environment. 

During each orbit at (3 = 0°, HALOE collects data for a maximum total of 15 minutes (7.5 min. during 
each sunrise and sunset), and is stowed during the remainder of that period. Assuming that only the volatile 
species accumulate on critical optical surfaces, the average mass flux rate <j> per orbit per species j may be 
estimated by 



In Eq. (2), linear behavior of <j>j is assumed between the two cases, and in Eq. (3) for conservative estimates, 
^>,aunaet = ^j.aunriaej even though satellite IM-area viewfactors are lower for sunset viewing periods. 

The probability that a contaminant molecule adheres to a critical surface is determined by its “sticking co- 
efficient,” a probabilistic parameter dependent on the combined effects of individual contaminant /surface ma- 
terial pairing, energy differences between contaminant molecules and the surface, and photopolymerization, 15 
among other things. Multiplying Eq. (3) by the number of orbits in the 35-month nominal HALOE lifespan 
and the area of the HALOE aperture, an upper bound for the total mass accumulated for each species may 
be calculated, assuming a sticking coefficient of unity. Neglecting the effects of highly-reactive O atomic 
flux on removing these deposits, such an assumption may be valid during data collection periods due to 
photopolymerization of volatile compounds on optical surfaces, but may considerably overestimate the mass 
accumulated while HALOE is stowed facing away from the Sun. Assuming values for the mass and volume 
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associated with each species, one can make a crude estimate for the depth of deposits at the aperture for 
each species. The total thickness becomes the summation of the contributions from each species. 

5.1»2 DSMC Resolution 

Figure 3 illustrates representative cross-sections of the inner domain cell network used for the present 
study after grid adaptation. Figure 3b shows the extent of the solar panel along with the rest of the vehicle. 
The body geometry is discretized and stored at the FCM level. Each element is characterized by a wetted 
surface area and a set of normal direction cosines. Figure 4 shows an XY-planar slice through the discretized 
geometry using a diagnosis utility developed for use with the DSMC code. In Fig. 4, pixels representing the 
surface are displayed along with outward direction cosines. Additionally, information on different levels of 
subcells surrounding the geometry can be represented. To simulate the UARS geometry, about 85,000 pixels 
have been used (pixel dimension % 5.0cm), and the inner domain computational grid is composed of 35,000 
cells. All outgassing and venting particles enter the flowfield from positions adjacent to the center of the 
pixel face closest to the outward normal of the outgassing pixel, about 0.2 percent the distance represented 
by a pixel length scale inside the flowfield volume. This has been done to eliminate position round-off error. 
These particles are given velocity components sampled from a drifting Maxwellian at the local outgassing 
surface temperature. 

5.2 FREESTREAM PARAMETERS 

Assumed ambient conditions at the satellite operating altitude of 600ifcm consisted of 100 percent atomic 
oxygen (0) at a nominal number density of 6xl0 12 molecules/m 3 . 16 At this altitude, the density may vary 
from this standard in either direction by roughly an order of magnitude due to solar activity. In addition, 
solar activity can cause appreciable helium concentrations to develop at this altitude. For conservative 
calculations, assuming 100 percent 0 maximizes its presence, and since 0 is highly corrosive to certain 
satellite materials, and may actually reduce the cumulative level of deposits on critical optical surfaces, 
instrument managers may be interested in its interactions with the satellite. Of course, ambient collision 
rates are affected by number density and species constituency, 17 however, using a single species to represent 
the freestream simplifies the simulation process, and it happens that UARS is operating in a period of low 
solar activity. These reasons may make the assumed level for rioo somewhat conservative. The assumed 
satellite orbital velocity V $ is 7500 m/sec. 
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5.3 QUTGASSING SPECIES. FLUX RATES. & TEMPERATURES 
5,3.1 Outgassing Species 

In preliminary calculations, 6 two ficticious outgassing species were fabricated for demonstration purposes. 
One, with a gram molecular weight (GMW) of 200, was emitted from all outgassing surfaces. The other, 
with GMW = 100, was emitted only from the Instrument Module (IM) vent, located above the CLAES 
experiment. 6 The GMW's are similar to those used in standard industry contamination codes, and are 
associated with heavy, volatile, organic compounds. 5 

Further investigation into likely types of outgassing materials covering UARS focused on four differ- 
ent items. These were: Mylar, used in the Multi-Layered Insulation (MLI) covering most of the satellite; 
Chemglaze Z306, a black paint covering portions of some instruments (as well as the interior of HALOE); 
S/13G/LO-V10, a white paint covering portions of some instruments as well as the back side of the SP; 
and adhesives possibly used to affix components to the front of the SP. A report investigating the types and 
mass fractions of volatile components outgassing from these materials at different temperatures identified 
compounds such as alkanes, alcohols, organic acids, and aldehydes. 18 Many components remained unidenti- 
fied; however, taking the highest mass-fraction species identified for the materials listed above, a ficticious 
molecule with GMW close to each calculated average has been generated. Moreover, a ficticious species 
associated with venting from internally-carried components has been created. The GMW assumed for this 
material is 100, following standard contamination procedures. 5 Values describing each of these species and 
their outgassing rates are listed in Table 10. 

Collision cross-sections for these species ratioed to that for monatomic oxygen have been assumed to scale 
with the 2/3 power of the associated mass ratio (see Table 10). 17 The mean viscosity-temperature power-law 
exponent u is assumed to be 0.75, a value associated with O, and no chemistry has been considered. This 
information is contained within the input gas species file used by the DSMC code, which is presented in 
Table 11. Due to lack of information, it is not clear how well the collision dynamics of heavy, long-chain, 
polyatomic organic molecules are modeled by these assumptions, but since temperatures are moderate and 
density levels are extremely low, collision rates are assumed to be negligible. 

All surfaces reflected particles diffusely, although evidence suggests incomplete momentum and energy 
accommodation better represent certain surfaces exposed to the vacuum of space for extended periods. 1 ^ 
Such behavior could affect the statistics of true fluxes encountered at the instrument aperture. 

5.3.2 Mass Flux Rates 


Measured mass flux rates associated with materials covering the instruments and subassemblies have been 
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obtained from Ref. 3, which lists data taken from a month-long bakeout test performed at NASA Goddard 
Space Flight Center on the completed UARS satellite at 100 °C. In that report, rates were assumed for 
units lacking measured data. It was noted that the rate associated with MLI covering the satellite bus 
varied from a high initial level to a lower, constant level by the end of the test. In this study the lower, 
constant value was assumed, since the HALOE aperture remained closed for the first month in orbit as a 
precautionary measure against the initial, transient levels of MLI outgassing. Unfortunately, no information 
was given on the relative location of instruments measuring these rates with respect to the satellite, and 
some subassemblies feature surfaces sporting more than one type of covering. Since only one rate was 
listed per unit, such information had to be supplemented through judicious interpretation of photographs 
depicting various UARS subcomponents to determine what surfaces were covered with which material. These 
photographs were supplied by HALOE instrument managers and personnel at Goddard Space Flight Center. 
The mass flux rates for surfaces and coatings deduced from Ref. 3 are presented in Table 10. 

In addition to outgassing surfaces, other gas-emitting sources have been modeled. The CLAES experiment 
is cryogenically cooled by gases contained in two separate reservoirs. 20-22 One reservoir contains Ne, the 
other C02. After cooling the instrument, these gases are expelled through vents. The CO 2 vent is located 
at the exit of a 1 in-diameter tube intersecting the CLAES shield plate, with direction cosines for the exit 
plane normal as (-.49, .64, .59) in the X, Y, Z coordinate system indicated in Fig. 1. The Ne vent is 
located at the end of a short, 2 in-diameter tube at the +Z - most section of the CLAES smaller cylinder, 
directed downward. The flow from these vent exits is assumed to be choked at an estimated temperature 
of T ex it — 207 K? 2 Mass flow rates for the CO 2 and Ne vents have been estimated at 3.65 and 8.39m<j/sec, 
respectively. In this study, these sources have been modeled using the jet-like vent plume routines for sonic 
orifices mentioned above. Calculated outgassing rates for the CLAES vents are listed in Table 10. 

Other vents, such as the IM vent, MLI blanket vents, and those associated with the MMS, released effluents 
at much lower flux rates, and have been modeled as regular surface-temperature dependent outgassing 
surfaces. Assumed GMW ' s and outgassing rates associated with these vents are listed in Table 10. 

5.3.3 Temperature Variation 


Mass flux rates given in Ref. 3 for the bakeout test have been adjusted to reflect levels associated with re- 
alistic operating temperatures through use of the following empirical relationship, formulated to characterize 
generic mass flux behavior of polymeric materials on Skylab 4,5 ( t in °C): 

M 


ft - 100 \ 

= exp HH 


<t>j(t = 100 °c) V 29 ; ^ 

This formula has been used to adjust all volatile organic compound outgassing rates. Although actual 
material outgassing rates probably deviate from Eq. (4), this formula is the only one available for use in 
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this investigation. Surface temperatures associated with various components have been estimated using data 
contained within Ref. 3, with one exception. UARS instrument managers provided solar panel operating 
temperatures which were measured in orbit. For /? = 0°, SP temperatures oscillate on orbit between -60°C 
and 30°C. The maximum SP temperature has been assumed for all cases, and values for other subassembly 
surfaces are listed in Table 7. 

5.3.4 Species Weighting Factors 

Upon commencement of the first DSMC computations, problems were encountered obtaining volatile fluxes 
at the HALOE aperture, because the outgassing rates associated with those species created flowfield number 
density levels orders of magnitude lower than those associated with the freestream and the CLAES vent 
species. 7 In a single timestep, about 3300 simulated monatomic oxygen molecules would enter the upstream 
boundary (area 25 110m 2 ), 30 simulated volatile, organic molecules of all types outgassed from the satellite 
(surface area 25 230m 2 ), and 200,000 particles exited the CLAES vents (total vent exit area = 2.5xl0 -3 m 2 ). 
The number flux of volatiles in the inner domain is therefore about two billionths of that for the CLAES 
vent species. Thus, by adjusting computational resources to adequately simulate the dominant C0 2 and Ne 
environment, keeping the computational size of the simulation large enough to simulate ~ 10 6 particles, the 
simulation would completely lose track of the volatile species whose development was the motivation for 
engaging in this investigation. 

The solution adopted has been to implement the species weighting factors discussed in Section 4.3. First, 
a simulation was run without the contributions of the CLAES vents in order to assess the ratio of volatile, 
organic species number density levels to the freestream. Another simulation included CLAES vent con- 
tributions in order to determine those density levels relative to those from the freestream as well. With 
information on the relative species number density levels, simulated fluxes of the various species have been 
altered using arbitrary SWF’s in order to obtain roughly equal orders of magnitude for numbers of each 
simulated species in the computational volume. The weighting factors used for O, Ne, C0 2 , and the volatile, 
organic species with GMW = 100, 150, and 200 are 40, 1, 1, 10 3 , 10<, and 5xl0 3 , respectively. 

5.4 CODE SETUP 

The methodology adopted for this study is to initially run the code for each set of conditions in a free 
molecule (collisionless) mode. Due to the large size of the satellite and the low thermal velocity of the heavy 
outgassing molecules, the approach to steady state is slow (N„ 25 10 4 At, where At » 7/isec). Figure 5 shows 
the number of particles in the inner and outer domains versus number of elapsed timesteps. Depending on 
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the type of problem being studied, the region where the curves flatten out may mark the onset of steady 
state. (While this rule-of-thumb may usually be adequate for measuring surface quantities, it is no panacea. 
For example, wakes behind bodies evolve over a long period of time without perceptible change in overall 
density levels.) Since exploratory calculations 6 ' 7 demonstrated that the global level of collisions was very 
low, the simulation is initially run in free molecule mode to reach an approximate steady state more quickly 
than if intermolecular collisions were allowed. After grid adaptation, each adapted cell contains roughly 
the same number of molecules (about 10 for the inner domain), as shown in Fig. 6. The code is then run 
with intermolecular collisions until a new steady state is reached. Subsequently, molecules incident on the 
HALOE aperture are counted, sampled, and written to a data file. 

For the inner domain, high spatial resolution is necessary, since UARS is very complex and quite large 
(~10m), and the HALOE aperture is only 22.9cm (9 m) in diameter. At the CCM level, the inner domain 
is composed of 74x46x104 equal subdivisions of 15.1cm on a side. Each FCM subcell consisted of pixels 
having lengths three times smaller than the CCM level, resulting in a surface resolution of 5.0cm. Since 
overall density levels are lower in the outer domain, the physical domain could be simulated with much 
coarser resolution. The outer domain is a cube measuring 21.2m on a side, composed of 70x70x70 subcells 
at the CCM level, with each FCM subcell having resolution equal to the CCM level. 

6.0 RESULTS AND DISCUSSION 


The wealth of information obtainable within a DSMC simulation will clearly be demonstrated in this 
section, as surface flux information is obtained along with flowfield behavior within the three-dimensional 
solution volume. First species number density contour maps will be presented, describing concentration levels 
surrounding the satellite, followed by results characterizing the fluxes incident on the HALOE aperture. The 
latter includes statistics concerning prior collisions, velocity distribution functions, and mass flux rates 
associated with simulated molecules impinging upon the aperture. This last item is manipulated to obtain 
a cumulative contamination deposit level. Finally, this cumulative level is compared to that of Ref. 3, and 
some code performance measurements are presented. 

6.1 SPECIES NUMBER DENSITY 

Figures 7-15 depict number density contour maps for each of the six gaseous species in the entire simulated 
flowfield surrounding UARS, ratioed to the freestream ambient number density, no, oo using a log 10 scale. 
This information is presented in two planes that capture information concerning the majority of the flowfield 
physics in the three-dimensional volume. One plane runs longitudinally, parallel to the main body axis, while 
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the other runs normal to it and the freestream direction. 

For the first set of figures concerning Case 1 (Figs. 7-12), a longitudinal solution plane has been chosen 
that runs through the vehicle coincident with its main body axis (X-axis), while the transverse solution 
plane coincides with the streamwise location of the center of the CLAES Ne vent in order to capture a 
cross-section of the plume emanating from it. These planes have been picked to present the overall physical 
environment. For the remainder of these number density maps, planes have been chosen that better reflect 
the environment in the vicinity of the HALOE instrument. 

The black outline surrounding UARS denotes the boundary separating inner and outer domains. While 
contributions of each species are presented separately, they co-exist simultaneously within the simulation. 
6.1.1 Case 1 

In Figs. 7a and 7b, contours for freestream species monatomic oxygen are presented for Case 1. In 
Fig. 7a, an expected buildup of O is observed on UARS windward side components. However, in initial 
DSMC calculations omitting the CLAES vents, 6 and in subsequent calculations prior to grid adaptation, 
when intermolecular collisions have been disallowed, the maximum density level peaked at about 25 times 
the freestream level on the cylindrical face of the MMS. 6 After grid adaptation and allowing intermolecular 
collisions to occur, the maximum density ratio there becomes roughly 50 times the freestream level. It seems 
apparent that this density level increase is due to interaction of the freestream with the gases emitted from 
the CLAES vents. Even though the global level of intermolecular collisions is low (about 130 collisions per 
timestep over 740,000 simulated molecules), the CLAES vent plumes create an effectively larger frontal area 
for UARS with respect to the freestream and scatter O out of the streamwise direction. An advantage of 
presenting these maps with a logarthmic scale is that the diminished levels of O in the wake of the vehicle 
are clearly evident. Figure 7b depicts the same set of contours viewed from the opposite side of the domain. 
The large, dark parallelogram jutting out from behind the transverse plane is the SP, canted 6° upward for 
this case. 

In Figs. 8 and 9, number density contour maps are depicted for CLAES Ne and C0 2 vent gases. In 
Fig. 8a, the plume developing from the Ne vent exit is clearly captured. In the transverse plane, one can 
observe a sizable amount of upward scattering although the vent is pointed directly downward. This is likely 
the result of diffuse surface scattering from surfaces of the CLAES experiment, due to its proximity to the Ne 
vent. Figure 8b reveals that very little of the Ne plume extends into the -Z regions. In Figs. 9a and 9b, the 
plume associated with the C0 2 vent is also clearly defined, although its axis is directed at oblique angles to 
both the streamwise and transverse planes. 

In Figs. 10-12, number density contour maps are presented for the outgassing species. In Figs. 10a and 10b, 
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the entire vehicle is surrounded by outgassing of species with GMW = 100- The streamwise plane contour 
map above UARS is dominated by outgassing from the IM vent (located just above CLAES on the +Z side 
of the satellite bus and angled 45° upward), although much of the concentration extending in front of the 
satellite is due to outgassing of adhesive material from the cell side of the SP (see Fig. 10b). The minor 
contribution to the flowfield from surfaces outgassing Chemglaze Z306 is shown in Figs. 11a and lib (GMW ; 
= 150). These clouds originate chiefly from surfaces located on NEPS and ZEPS, although the interior of 
HALOE outgasses this material through its aperture. In Figs. 12a and 12b, GMW = 200 concentrations are 
presented. In the streamwise plane, the highest density levels are found near the HGA dish, although most 
of the outgassed material originates from the back side of the SP (Fig. 12b). 

Results for this identical case are shown using two other planes that better depict the HALOE environment. 
The HALOE instrument is located in the center of each figure. In Fig. 13a, a buildup of O is clearly seen in 
front of HALOE and its mount. One also notices that the highest levels are found before the SP as it plows 
through the rarefied atmosphere. In Fig. 13b, even though the vast majority of Ne is directed downward on 
UARS’ other side, there is a slight presence of that species in the vicinity of HALOE through a combination 
of surface and intermodular collisions. As for CO 2 , in comparing streamwise plane results for Fig. 13c 
with Fig. 9b, where the former is displaced approximately lm closer to HALOE than the latter in Fig. 9b, 
density levels for that species have abated by roughly an order of magnitude. No CO 2 flux was detected at 
the HALOE aperture. 

Comparing streamwise planes from Fig. 13d with Fig. 10b for the GMW = 100 material, it becomes 
clear that the influence of the IM vent is greatly diminished over that lm separation distance, and the local 
flowfield is dominated by the influence of the SP. In Fig. 13e, the contribution of the GMW = 150 material 
is globally quite minor, yet there is a relatively high concentration of it at the HALOE aperture, thanks 
to its own contribution. Finally, the concentrations of GMW — 200 material at the aperture are shown in 
Fig. 13f to be about an order of magnitude higher than for the previous figure, which is a consequence of the 
high concentrations of this material from the back side of the solar panel and the — Z side of the Instrument 
Module. 

6.1.2 Case 2 

Planes chosen to represent the Case 2 flowfield solution (Figs. 14a-14f) are identical to those depicted in 
the last subsection, emphasizing the HALOE environment. The number density contour maps for this case 
are similar to those presented in the preceeding two paragraphs. Besides the obvious repositioning of the SP 
and the HALOE instrument (-23° inclination versus 6° above), the most striking difference is in the map for 
the GMW = 200 material (Fig. 14f vs. Fig. 13f). The difference in SP inclination accounts for a dramatic 
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decrease of this material in the vicinity of HALOE. 

6.1.3 Case 3 

The general features presented in Figs. 15a through 15f are similar to the other two cases. Although 
the SP front face and HALOE no longer face the freestream, the SP rear face now opposes it. In Fig. 15a, 
the solar panel still impedes the motion of freestream 0, but the position of the transverse planar solution 
cut reveals the wake behind it. In Fig. 15b, Ne is still present around HALOE, but HALOED orientation 
25° away from the body axis toward -Z effectively shields its aperture from this gas. No Ne or C0 2 was 
detected at the aperture for this case (Figs. 15b and 15c). 

Because the SP is roughly aligned with the +X direction, the GMW = 100 contour density map in 
Fig. 15d shows a striking difference from the analogous Case 1 solution depicted in Fig. 13d above, where 
it is roughly aligned in the opposite direction. There is a reversal of which portions of the flowfield are 
densely and sparsely populated with this material. Even though HALOE is stowed, its viewfactor of the SP 
approaches that for Case 1 (Dsp|ca«e 3 = 0.412 steradians versus f2sp|cas« l = 0.527 steradians). It will be 
shown shortly that in this configuration, there is still a surprisingly significant level of contaminant flux at 
the instrument aperture, mostly due to the SP. Fig. 15e shows the contribution from GMW = 150 material 
for this case, with local concentration maxima emanating from the HALOE aperture, NEPS, and ZEPS. 
Finally, in Fig. 15f, number density contours for the GMW = 200 material are presented, where once again, 
a reversal in high and low concentrations of this material are evident compared to Case 1 results in Fig. 13f. 

6.2 CAPTURED MOLECULES-PRIOR COLLISIONS 

There are a variety of statistics one can glean from simulated molecules striking the HALOE aperture. 
Two sets of information are presented in this section and the next. First, for each case, the history of each 
impinging simulated molecule is presented with respect to number of prior collisions it has suffered before 
reaching the aperture. Second, histograms are generated, representing the AT, Y } and Z velocity probability 
distribution functions (PDF) for these same particles. 

6.2.1 Case 1 

In Figs. 16a- 16c, plots are presented depicting the number of surface collisions prior to reaching the 
aperture for each species in each case. For Case 1, depicted in Fig. 16a, nearly 90 percent of O atoms reach 
the aperture directly from the freestream with no prior surface collisions. In contrast, all Ne atoms suffer at 
least two surface collisions before reaching this location. Outgassing molecules with GMW =100 and 150 
tend to encounter multiple surfaces before arrival, while most GMW = 200 material deposits directly from 
the back side of the SP and the -Z side of the IM. 
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6.2.2 Case 2 

In Fig. 16b, Case 2 collision history statistics for O are similar to those for Case 1. However, Ne atoms 
collected encountered at least five prior surface collisions before impinging on HALOE. This could be because 
the solar panel lies behind it in this run, and in its visual absence (zero viewfactor from the HALOE 
perspective), many other surfaces are required to facilitate the presence of neon. Again, outgassing molecules 
tend to encounter multiple reflections before arriving at the aperture, except for the GMW = 200 material. 
This material tends to collect from direct viewing of the — Z side of the IM. 

6.2.3 Case 3 

Figure 16c shows that no O atoms directly encounter the aperture; most reflect once off the SP instead. 
Also, due to HALOE instrument orientation, no Ne is detected at all. GMW = 100 and 150 materials also 
tend to reflect once off the SP before impinging, but as in Case 1, GMW = 200 material emitted from the 
back side of the SP tends to arrive directly. 

These figures do not reveal all information pertinent to O , whose distributions contain a rather long tail 
of multiply-reflected particles. The highest number recorded is 1 14 prior surface collisions for Case 1 . 

6.3 CAPTURED MOLECULES- 

VELOCITY PROBABILITY DISTRIBUTION FUNCTIONS 
6.3.1 Case 1 

Figure 17 depicts Case 1 histograms representing the X , 7, and Z PDF for particles striking the aperture. 
Only Ne contributes any return flux; all other species have been deposited as the result of direct flux. Be- 
ginning with the X , or streamwise direction in Fig. 17a, notice that PDF(V X ) for O is bimodal. One mode, 
a full Maxwellian distribution centered at V x = V $ , represents molecules having suffered no previous surface 
collisions. The other mode, corresponding to a drifting Maxwellian distribution near V x = 0, represents 
molecules having previously collided with the satellite at least once. Figure 17b shows PDF behavior for 
V T near the origin for better detail. In this figure, the volatile outgassing species have a streamwise PDF 
corresponding to drifting Maxwellian distributions whose shapes differ from one another due to individ- 
ual GMW's. These distributions have been checked and compare favorably with theoretical calculations, 
but for clarity’s sake the theoretical functions for PDF(1/ X ) have been omitted from the figure. However, 
the Ne distribution does not conform to any simple theoretical model because the majority of Ne samples 
had experienced intermolecular collisions. In Figs. 17c and 17d, probability density functions are presented 
for the 7 and Z, or transverse, directions. The O distributions behave as expected, as they are dominated 
by contributions from the freest ream mode, and compare quite well with theory. Due to the complicated 
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nature of the geometric setup, no comparisons between these histograms and theoretical PDF’s have been 
attempted in these directions. 

6.3.2 Case 2 

Figure 18 depicts Case 2 velocity PDF histograms. The X direction plot (Fig. 18a), generally appears 
similar to that for Case 1. However, the origin detail displayed in Fig. 18b shows that some portion of the 
flux has V x < 0. Figs. 18c and 18d show that although O behaves as in Case 1, other species PDF’s tend 
to be narrower and more peaked than in the previous case. This may be because overall viewfactors for 
HALOE are lower than Case 1, and impinging molecules come from a smaller, more select group of surfaces. 

6.3.3 Case 3 

In Case 3, there is no contribution from directly-impinging, freestream O atoms, so the distinct, bell- 
shaped Maxwellian mode associated with that species for Figs. 17a and 18a is not present in Fig. 19a. 
Although the mass flux associated with O still dominates the HALOE aperture environment, it is nearly 
two orders of magnitude lower than for the first two cases. The contribution from GMW = 150 molecules 
come entirely from the aperture itself, where emitted molecules reflect off the SP back towards HALOE. 
In Figs. 19b and 19c, most of the O flux has < 0 and V z > 0, indicating that it is constituted mostly 
from simulated atoms scattered downward and toward UARS from the SP, since that is the only assembly 
occupying -Z space relative to HALOE. 

6,4 MASS FLUX RATES 

Individual species mass flux rates impinging on the HALOE aperture are presented in Table 12. Accuracy 
of these results rests heavily on the reliability of values for surface outgassing flux rates at the bakeout 
temperature, the use of Eq. (4) to adjust those rates for surface temperature, and the relative orientation 
and proximity of HALOE to the contributing surfaces. The rate for O has been subdivided into those atoms 
coming directly from the freestream and those having suffered prior surface collisions. Since a substantial 
portion of impinging Ne is attributed to return flux, while the rest reaches HALOE through surface collisions 
alone, two entries have been listed for Ne as well. Notice that fluxes from the outgassing species are orders of 
magnitude less than the others. Without using SWF’s, much longer run times would be required to achieve 
a reasonable sample sizes — possibly by a factor of many thousand. 

Tables 13a-13c relate the relative contributions of different surfaces to the mass flux rates presented in 
Table 12. The SP is the chief progenitor of the HALOE contamination environment for all three cases. This 
is because it apparently attains the highest temperatures of any surface expected on the satellite for = 0° 
(see Table 7), it is large, and it is located relatively near the HALOE aperture. The -Z side of the IM 
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(second largest contributor), also represents a large outgassing surface near HALOE, but its contributions 
are diminished because its outgassing rates are lower than those of the SP, and its surface normal does not 
oppose the HALOE aperture. 

Flux rates for Case 2 are significantly less than Case 1 due to reduced viewfactors of the main contributing 
outgassing surfaces. Included in Tables 13a-13c are viewfactors ([solid angle Q] = steradians out of 2ir) of 
the SP and the — Z side of the IM experienced by the HALOE aperture. Ironically, Case 3 demonstrates 
that while a stowed HALOE instrument is exposed to very little direct flux from the bulk of UARS, it is 
sometimes prone to relatively high viewfactors of the continuously-rotating SP, periodically enhancing the 
pernicious contamination environment surrounding HALOE. 

Also enhancing this environment is a contribution from the Chemglaze-painted surfaces within its interior, 
whose effluents are allowed to escape out through the aperture. Some of these molecules strike other surfaces 
and return to the HALOE aperture, and are counted along with the other contaminants. 

6.5 CUMULATIVE DEPOSIT LEVELS 

Eqs. (2) and (3) may be used to estimate an upper-bound cumulative total level of deposits at the HALOE 
aperture. Assuming that only volatile, organic molecules accumulate (with a sticking coefficient of unity), 
and each occupies the space of a cube of length /, where / equals the particular species diameter listed in 
Table 10, a deposit depth of « 3550A is estimated. If standard contamination procedures are followed, and / 
= 10A is assumed regardless of species diameter, 3 that estimate increases to ss 13,350A. If the contribution 
HALOE receives while stowed is disregarded (sticking coefficient of zero), these estimates fall to ss 650A and 
« 2650A, respectively. 

Tables 14 and 15 show the species-level contributions from each run to the final, averaged estimates quoted 
above. It is interesting to observe how GMW factors into the second estimate described, which is a standard 
procedure in the contamination community. 3 A greater number of lighter GMW species molecules are required 
for a mass flux equal to that of a heavier GMW species. Since the deposit depth of contaminant molecules 
only depends on the number that accumulate, it becomes important to accurately assess the GMW's of 
the accumulated molecules. Also, the true sticking coefficient is not known, and can be much less than 
unity under certain conditions (small gas/surface temperature differences, absence of direct sunlight, etc.). 15 
Finally, the effects of scrubbing from the high monatomic oxygen flux on reducing the level of accumulants 
is not clear. 

The effect of this predicted cumulative buildup on instrument performance is not clear. At this time, 
UARS has been in orbit for over two years, and the instrument has shown negligible loss in performance. 
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HALOE instrument managers have provided figures quantifying performance degradation in terms of signal 
strength attenuation for different species concentration measurements over the first year of operation. Signal 
voltages obtained over the first two months of operation were used to compute a standard expected signal 
voltage strength per channel, with instrument temperature corrections. Subsequent voltage strengths were 
measured and compared to the standard. Figures 20a and 20b show signal loss correlations for HF and HCl 
channel measurements, respectively. For initial standard voltages of » 3.3 volts and ss 3.45 volts, after one 
year’s operation, HF and HCl signal strengths have apparently been attenuated by 0.02 volts and 0.004 volts, 
respectively. This level of voltage attenuation is indicative of measurements made for the other channels 
used by HALOE as well. This may suggest that the concentrations of contaminants in the vicinity of the 
HALOE instrument have been overestimated in this as well as in previous studies, or that the effects of 
contamination from the types of outgassing materials carried on board UARS are relatively benign. Perhaps 
better understanding of volatile GMW's, outgassing rates, temperature dependence on outgassing rates, 
gas/surface sticking coefficients, the cleansing effects of atomic oxygen, and light attenuation and scattering 
due to outgassing contaminants will come from the current examination of samples carried aboard the Long 
Duration Exposure Facility (LDEF). 23 

6.6 COMPARISON WITH PREVIOUS ANALYSES 

It is difficult to make detailed comparisons with the results described in Ref. 3 for a number of reasons. 
First, in Ref. 3, the orientations of the SP and HALOE were decoupled. The SP cellside surface was pointed 
downward at HALOE, minimizing the distance between the two objects. Meanwhile, HALOE was allowed 
to assume various orientations associated with data collection periods as well as possible stow positions. In 
reality, HALOE would be stowed pointing away from the SP when its cell side faces downward. Second, 
the earlier contamination analyses 3 assumed a constant solar array surface temperature of 100 °C. This 
investigation benefitted from on-orbit data measured directly on the SP which indicate that 0 = 0° surface 
temperatures are cyclic with each orbit over a range of -60 *C to 30 °C. For the 30 °C value used in this 
study, Eq. (4) predicts that the SP would outgas only 9 percent of the value expected at 100 °C. Third, 
the earlier analyses assumed only one outgassing species ( GMW = 100) to relate cumulative depth to mass 
fluxes. Obscuring detailed comparisons further, the earlier analyses could not take into account the fact 
that after the Goddard bakeout test, MLI blankets covering the instrument module surface on the — Z side 
had themselves been covered with panels of S/13G/LO-V10 (which outgasses at mass flux rates one to two 
orders of magnitude higher than MLI see Table 10). This had the impact of increasing assumed flux levels 
from those surfaces to HALOE. Finally, the geometry used in Ref. 3 was somewhat less detailed than in the 
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present study. Taking these differences into consideration, it seems that the results of this investigation tend 
to confirm those contained within Ref. 3 nevertheless. 

6.7 CODE PERFORMANCE 

The DSMC code is run simultaneously on two Sun SPARC 2 workstations with 64MBytes memory (RAM) 
each. Statistics concerning memory requirements, number of particles, weighting factor FNUM, timestep At, 
and overall computational volumes for the cases studied are listed in Table 16. 

Regarding values for CPU time spent processing each region, it should be noted that FNUM and At could 
have been adjusted to create better load balancing between the inner and outer domain processors. This 
was not a major goal of the study, however. 

Computational performance for simulations using particle methods is often measured in terms of average 
CPU time taken to process the state of an average particle per timestep. Using this metric, levels of 
approximately 29 /jsec/particle/timestep and 44 /jsec/particle/timestep have been obtained in the outer and 
inner domains, respectively, with the difference chiefly due to handling outgassing surfaces in the latter. 
This is a fairly exotic application for DSMC however, and it may be difficult to compare this performance 
estimate to other results in a straightforward, detailed manner. 

7.0 CONCLUDING REMARKS 

This report describes the development and application of a DSMC code to study contamination problems 
on satellites and estimate the effect of outgassing on sensitive instruments. Discussion is included on the 
development of a number of specially-developed preprocessor routines to incorporate detailed information 
regarding geometry, outgassing rates and temperature data, different types of outgassing sources, and the use 
of species weighting factors to enhance the statistical representation of important trace species. In estimating 
the amount of volatile organic molecular deposits collected at the HALOE aperture, special attention has 
been given to accurately account for outgassing species characteristics, surface temperatures, and geometric 
orientation of UAES subassemblies deemed critical to this endeavor. Using the best information available, 
and neglecting the scouring effects of monatomic oxygen, an upper bound on the cumulative contaminant 
deposition level at the HALOE aperture has been estimated at 13,350A. If there is negligible tendency for 
volatile species to accumulate while HALOE is stored, this value falls to 2650A. Notwithstanding consider- 
ation of important differences in geometric orientation, outgassing species weights and flux rates, and other 
specific influential input data, it appears that these results tend to confirm the earlier work described in 
Ref. 3. 
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Although the model for outgassing contained in DSMC has often been compared to other methods, 910 ' 24 it 
could be instructive to make direct comparisons between results obtained using the DSMC method described 
in this report for a simple geometry over different ambient conditions with results from industry-standard 
codes and experimental data when they become available. Also, this study would benefit from better 
characterization of outgassing materials in terms of average molecular weights, mass flux rates, temperature 
dependence on mass flux rates, etc. 


REFERENCES 


1 Gordley, Larry L., Russell, James M. Ill, and Hesketh, W. Donald, “Calibration, Testing and Orbital 
Performance of the Halogen Occultation Experiment (HALOE) on the Upper Atmosphere Research Satel- 
lite,” AIAA Paper 93-0128 , Jan. 1993. 

2 Robertson, S. J., “BGK Model Solution of Back Scattering of Outgas Flow from Spherical Spacecraft,” 
Progress in Astronautics and Aeronautics, Vol. 51, Part I, July 1976. 

3 Hatef, J. B., “Predicted Contamination of the HALOE,” Goddard Space Flight Center Transmittal 
Memo 7 32-91-011, May 1991. 

Anon., Shuttle/Payload Contamination Evaluation Program— The Space Computer Program User’s 
Manual,” Martin Marietta Aerospace, MCR-81-509 , Jan. 1981. 

8 Anon., “MOLFLUX-Molecular Flux User’s Manual,” Structures and Mechanics Division, Materials 
Branch, Rev. 1, NASA Lyndon B. Johnson Space Center, Houston, Texas, JSC-22496, Feb. 1989. 

6 Rault, D. F. G., and Woronowicz, M. S., “A Study of Instrument Contamination Aboard UARS Using 
Direct Simulation Monte Carlo,” AIAA Paper 93-0724, Jan. 1993. 

7 Woronowicz, M. S., and Rault, D. F. G., “On Predicting Contamination Levels of HALOE Optics 
Aboard UARS Using Direct Simulation Monte Carlo,” AIAA Paper 93-2869, July 1993. 

8 Bird, G. A., Molecular Gas Dynamics and the Direct Simulation of Gas Flows, 1— ed., Clarendon Press, 
Oxford, 1994. 

Rantanen, R. O., and Gordon, T. D., “Contaminant buildup on ram facing spacecraft surfaces,” Pro- 
ceedings of SPIE, Optical System Contamination: Effects, Measurement, Control, ed. Peter Glassford, 

Vol. 777, May 1987. 

Justiz, C., Sega, R., Dalton, C., and Ignatiev, A., “Return Flux Contamination of an Outgassing Space- 
craft in Low Earth Orbit,” AIAA Paper 93-0725, Jan. 1993. 



NASA TM 109069 25 


11 Bird, G. A., “Application of the Direct Simulation Monte Carlo Method to the Full Shuttle Geometry,” 
AIAA Paper 90-1962 , June 1990. 

12 Rault, D. F. G., “Aerodynamics of Shuttle Orbiter at High Altitudes,” AIAA Paper 93-2815 , 

July 1993. 

13 Rault, D. F. G., Wilmoth, R. G., and Bird, G. A., “An Efficient DSMC Algorithm Applied to a Delta 
Wing,” AIAA Paper 91-1316 , June 1991. 

14 Beguelin, Adam, Dongarra, Jack, Geist, Al, Manchek, Robert, and Sunderam, Vaidy, “A Users’ Guide 
to PVM Parallel Virtual Machine,” ORNL/TM-1 1826 , Oak Ridge National Laboratories, July 1991. 

15 Bettini, Ronald G., “Shuttle environment effects on coated mirrors,” Proceedings of SPIE, Optical Sys- 
tem Contamination: Effects, Measurement, Control , ed. Peter Glassford, Vol. 777, May 1987. 

16 U. S. Standard Atmosphere , 1976, U. S. Government Printing Office, Washington 25, D.C., 1976. 

17 Vincenti, Walter G., and Kruger, Charles H., Jr., Introduction to Physical Gas Dynamics , 1— ed., 
Robert E. Krieger Publishing Company, Malabar, Florida, 1982. 

18 Glassford, A. P. M., and Garrett, J. W., WRDC-TR-89-4 1 14 1 “Characterization of Contamination Gen- 
eration Characteristics of Satellite Materials,” Materials Laboratory, Wright Research and Development 
Center, Air Force Systems Command, Wright-Patterson Air Force Base, Ohio, Nov. 1989. 

19 Fredo, Richard M., and Kaplan, Marshall H., “Procedure for Obtaining Aerodynamic Properties of 
Spacecraft,” Journal of Spacecraft and Rockets, Vol. 18, No. 4, 1981, pp. 366-373. 

20 Bortner, Dr. M. H., “Effect of Vented CO 2 and Neon on Performance of UARS Instruments,” General 
Electric Astro Space Division, Valley Forge Space Center Progam Information Release, U-l K21-UARS- 
853 , June 1987. 

21 Thames, M., “Spatial Density Distributions for UARS CLAES Ne and C O 2 Vents,” General Elec- 
tric Astro Space Division, Valley Forge Space Center Progam Information Release, U-1K20-U ARS-878 , 
July 1987. 

22 Neste, S. L., “Effect of CLAES Venting on UARS Performance,” General Electric Astro Space Division, 
Valley Forge Space Center Progam Information Release, U-1K21-U ARS-919 , Aug. 1987. 

23 “LDEF-69 Months in Space, Second Post- Retrieval Symposium,” sponsored by LDEF Science Office, 
NASA Langley Research Center, and the American Institute of Aeronautics and Astronautics, ed. Arlene 
S. Levine, NASA CP 3194 , June 1992. 

24 Fan, C., Gee, C., and Fong, M. C., “Monte Carlo Simulation of Molecular Flux on Simple Spacecraft 
Surfaces due to Self- and Ambient-Scatter of Outgassing Molecules,” AIAA Paper 93-2867 , July 1993. 



Geometry Definition for X Program: 


26 NASA TM 109069 




dir. cosine 
radius at A 
radius at B 



Geometry Definition for X Program: 

Cylinders 


NASA TM 109069 27 



09 

3 

cd 

X 

£ 


< n 

•S ^ 

o o 

o< a, 


OQ 


cG 


•5 *3 

o cG 


a In T3 


0) 

a 

• rH 

O 43 

t a 

a> 


< PQ 


< PQ 


a> 

c 

• 

CO 

o 

a 


a> 

a 

> ^ 

CO 

o 

o 


c fla 

•l-H -l-H td l_i 

O O CG .B 

ft a Stjtj 


a 

o 


fl 

IB 

© 

Q 


0) 

X* 

c 


>> 

u 



co 


T3 



0Q 


T3 



a) 


£ 


rS 

OQ 



IS 

■5 

s> 

0) 

0) 

o> 

0) 

a 

C8 

s? 

T3 


T3 

.g 

^g 

g 





o 

u 

o 



Geometry Definition for X Program: 

Plates 


28 NASA TM 109069 






Geometry Definition for X Program: 

Rings 


NASA TM 109069 29 



GO 

3 

CO 

•p* 

u 

> 


u Ph 

co co 

C ^ ^ 

C *H * t— < 

•atf-o 

2 ts « 

& * * 


<N 


<N 


0) 

c 

CO 

o 


0) *H *H 

0) 0J 

0 ) 3 
o .H O T3 


a 

’3 aj 

ft ^ 

U4 QQ 

5 .2 

0 
o 


a> 

fi 

• H 

CO 

o 

a 


2 & 


0) -r* 

43 £ 

13 13 

43 43 


a 

o 


a 

tc 


4) 


fi 


bo 

o.g 

•c ~ 

-M U 

c « 

a> *73 
cj » 

3 g 

o .H 

o o 


“1 ft 

•c|| 

h § So 

J 2 o c 

3*5 $s 

•S '% 2 2 


S 

ee 

fc 


t 

a 


BS 


ffl 

be 

c 


• H 

Ph 



Geometry Definition for X Program: 


30 NASA TM 109069 
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Table 2: Geometry Input Data Spreadsheet 



Data 

Description 

27 347 


no. of data planes, quadrics 

0.0254 


scale 

100. 1. 

-20. 

origin of axis 

-80. 327. 

-100. 110. -515. 90. 

parallelpiped around body 

100. 1. 

-20. 

center of gravity 

10000. 300. 


reference area and length 

‘Main block top plate ABCD1’ 

Quadric 1 

1 


component to which this quadric belongs 

PlateB 



1 


flow is outward, inward, or on both sides 

3. 48.1 

-42.5 

pointA 

218. 48.1 

-42.5 

pointB 

3. 48.1 

2.1 

pointC 

218. 48.1 

2.1 

pointD 

100. 50. 

-20. 

exterior point 

0. 0. 0 

0 . 

displacement of limit planes (in pixels) 

‘CLAES small cylinder’ 

Quadric 15 

3 


component to which this quadric belongs 

CylinderA 



1 


flow is outward, inward, or on both sides 

64.5 2.7 

30.6 

pointA 

97.1 2.7 

30.6 

pointB 

20.0 


radius 

0. 0. 


displacement of limit planes (in pixels) 

0. 65. 

-30. 30. 0. 60. 

F3I limits 

‘CLAES Ring’ 

Plane 1 

3 


component to which this quadric belongs 

RingA 



64.5 2.7 

30.6 

pointA 

1. 0. 

0 . 

direction cosine perpendicular to plane 

20.0 


radius A 

25.9 


radius B 

‘CLAES small cap’ 

Quadric 16 

3 


component to which this quadric belongs 

SphereE 



1 


flow is outward, inward, or on both sides 

97.1 2.7 

30.6 

pointA 

-1. 0. 

0. 

direction cosine perpendicular to plane 

20.0 


radius 

0.5 


displacement of limit plane (in pixels) 

64. 118. 

-20. 25. 0. 55. 

F3I limits 
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Table 3: Input Data Listing — Geometric Model (29 pgs.) 
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Table 4: Input Data Listing — Subassembly Rotation Data 
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Table 5: Typical Elemental Outgassing Data 


Data 

Description 

‘Main block top plate ABCDT 

Quadric 1 

1 

number of outgassing surfaces on object 

PlateB 

shape of outgassing surface 

1 0 

flow direction code, outgassing index, 

(0 if entire surface has one rate, or no. of defining points) 

0. 1. 0.6264E+17 100. 

mole fractions of outgassing species, 

total outgassing flux [molecules/m 2 /sec], ref. temp. [°C] 

‘Main block side Zmax CDGHT 

Quadric 2 

2 

number of outgassing surfaces on object 

PlateB 

shape of outgassing surface #1, entire surface 

1 0 

flow direction code, outgassing index, 

(0 if entire surface has one rate, or no. of defining points) 

0. 1. 0.6264E+17 100. 

mole fractions of outgassing species, 

total outgassing flux [molecules/m 2 /sec], ref. temp. [°C] 

PlateB 

shape of outgassing surface #2, vent above CLAES 

1 4 

flow direction code, outgassing index, 
(number of defining points) 

1. 0. 0.7649E+20 100. 

mole fractions of outgassing species, 

total outgassing flux [molecules/m 2 /sec], ref. temp. [°C] 

12. 43.1 2.1 

pointA for vent 

60. 43.1 2.1 

pointB for vent 

12. 45.1 2.1 

pointC for vent 

60. 45.1 2.1 

pointD for vent 
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Table 6: Input Data Listing — Outgassing Information (19 pgs.) 
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Disk 

' HRDI telescope mount knob' 1 o 

1 Number of different outgassing surfaces on object 0. 0. 0. 1. 0. 0. 4.17E+14 100. outgas rate 

Disk <Number/m2/sec) , ref temp(C) 
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Number of different outgassing surfaces on object 'SSPP box ABEF14 ' 
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(Number /m2 /sec) , ref temp(C) 



Number of different outgassing surfaces on object {Number/m2/sec} , ref t«np(C) 
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of different outgassing surfaces on object 'WINDII base, GEM18' 



of different outgassing surfaces on object 'PEH-AXIS base, CDGH19 ' 
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behind antenna K0LP21' 


74 NASA TM 109069 




NASA TM 109069 75 




76 NASA TM 109069 




(Number /m2 /sec) , ref temp(C) 
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Number of different outgaesing surfaces on object (Number/m2/sec) , ref tanp(C) 




0. 0. 0. 1. 0. 0. 3.62E+14 100. outgas rate 'HALOE box DFJL2 8 ' 

(Number/m2 /sec I , ref tamp(C) 0 Number of different outgoing surfaces on object 

'HALOE mount V27' 'HALOE box MNOP28' 

1 Number of different outgassing surfaces on object 1 Number of different outgassing surfaces on object 

RingA PlateB 
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0. 0. 0. 0. 1. 0. 3.09E+15 100- outgas rate 

(Number /m2 /sec), ref temp(C) 'HALOE mount ACFG4 
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Number of different outgassing surfaces on object 
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Table 7: Input Data Listing — Surface Temperatures (3 pgs.) 
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Table 8: Typical Elemental Sonic Orifice Data 


Data 


Description 

0.0254 


scale 

0 . - 1 . 

0 . 

jet direction cosine, Ne vent 

71.0 1.7 

53.0 

jet center 

1 . 


jet radius 

3.28E+20 


jet density [molecules/m 3 ] 

377. 


jet velocity [m/sec] 

207. 


jet exit temperature [A] 

o 

o 

o 

0 . 0 . 

species fraction in jet 

-0.4924 0.6428 

0.5868 

jet direction cosine, CO 2 vent 

65.1 36.6 

26.2 

jet center 

0.5 


jet radius 

4.04E+20 


jet density [molecules/m 3 ] 

244. 


jet velocity [m/sec] 

207. 


jet exit temperature [A'] 

o 

o 

o 

0 . 0 . 

species fraction in jet 




data to model gas around UARS 
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Table 9: Input Data Listing — General Preprocessor Input 
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Table 10: Out gas sing Species Information 


Species 

GMW 

d 

[A] 

Mass Flux Rates 
at lOOT? [gm/cm^/sec] 

ficticious organic volatile species 
representing products from MLI 

100 

5.33 

4.33x 10 - ' 3 
to 6.01 x 10 -12 

ficticious organic volatile species 
representing products from solar array adhesives 

100 

5.33 

1.83x 10 -10 

ficticious organic volatile species 
representing products from IM and MLI vents 

100 

5.33 

1.60x 10 -7 (IM) 
1.22x10-® (MLI) 

ficticious organic volatile species 
representing products from MMS vents 

100 

5.33 

3.67xl0 -11 
to 1.22x10-'° 

ficticious organic volatile species 
representing products from Chemglaze Z306 

150 

6.33 

7.70x10-" 

ficticious organic volatile species 
representing products from S/13G/LO-V10 

200 

6.96 

4.30x10-" 
to 5.39x 10 - " 

neon vented from 
CLAES experiment 

20.18 

2.72 

4.14x 10 -4 * 

carbon dioxide vented from 
CLAES experiment 

44.01 

5.41 

7.20x 10~ 4 * 


* rate at operating temperature T ex j t = 207R 
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Table 11: Gas Species Data 


0 

7. 

30000 

0. 5 
0 

6 

2 

132 

7500. 

0. 

0. 

6 . e+12 1000. 

1 . 0 . 0 . 0 . 0 . 0 . 

300. 

300. 

2 2 2 0 2 2 
1 . 1 . 1000 . 

20 . 10 . 10 . 10 . 10 . 10 . 
350000 
2 

3. 0.28 

1 . 

2 50 

2 


spatial resolution switch: body-0, gas-1 
max no. of Megabytes for body defn. 
ratio to compute mean no. of subcells/cell 
ratio of cells near body 
initially uniform freestream or vacuum 
total no. of molecular species 

no. of jet species (subset of total no. of species) 

GAS?. FOR file ( species -dependent gas info.) 

free stream velocity 

angle of incidence alpha (deg) 

roll angle phi (deg) 

freestream density, temp. 

fraction of each species in freestream 

surface temp. 

mean expected temp, in flow 
boundary conds . at -X, -Y, -Z, +X, +Y, + Z 
safety factors for max no. of molecules, 
timestep, FNUM 

clearances for -X, +X, -Y, +Y, -Z, +Z edges 
no. of subcells within computational domain 
scale ratio between inner & outer domains 
anticipated density enhancements: peak, mean 

(used to estimate cell and subcell sizes) 
fraction of total volume to be used 
no. of timesteps between sampling, 

no. of samples between saved states 
parameter used to size pixel element arrays 


0 - 2 - - 
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Table 12: Mass Flux Rates Intercepted at HALOE Aperture 


Species 

Mass 

Case #1 

Flux Rates [gm/cm : 
Case #2 

! /sec] 

Case #3 

monatomic oxygen from 
freestream distribution 

1.2xl0 _I ° 

l.OxlO" 10 

< 4.2x 10 -14 * 

monatomic oxygen from 
surface-accommodated distribution 

9.3xl0~ 12 

7.5x 10~ 12 

5.4xl0 -12 

neon vented from CLAES 
experiment (return flux) 

1.9xl0 _u 

5.9xl0 -12 

< 2.1 x 10~ 12 * 

neon vented from CLAES 
experiment (direct flux) 

1.3x10-" 

5.9xl0 -12 

< 2.1xl0 -12 * 

carbon dioxide vented from 
CLAES experiment 

< 5.8xl0~ 12 * 

< 6.5xl0“ 12 * 

< 4.6x 10 -12 * 

ficticious organic volatile species 
representing all GMW = 100 materials 

9.2xl0 -14 

1.8xl0 -13 

8.4x 10 -14 

ficticious organic volatile species 
representing GMW =150 material 

3.9x 10 -15 

l.lxlO -14 

1 .4x 10 -14 

ficticious organic volatile species 
representing GMW = 200 material 

* i 

4.3xl0 -13 

1.8xl0 -14 

2.7xl0 -13 


* no samples collected 


Table 13a: Subassembly Contributions to 

Cumulative Contaminant Deposition (Case #1) 


Subassembly 

Percent Coni 
GMW = 100 

tribution to Conta 
GMW = 150 

minant Flux 
GMW = 200 

Solar Panel 
(ft = 0.527 steradians) 

56 

— 

81 

Instrument Module 
(ft = 0.719 steradians) 

13 

— 

19 

HALOE 

25 

73 

— 

NEPS j 

3 

27 

— 

ISA MS 

3 

— 

— 
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Table 13b: Subassembly Contributions to 

Cumulative Contaminant Deposition (Case #2) 


Subassembly 

Percent Con 
GMW = 100 

tribution to Conta 
GMW = 150 

minant Flux 
GMW = 200 

Solar Panel 
(ft = 0 steradians) 

92 

— 

— 

Instrument Module 
(ft = 0.292 steradians) 

— 

— 

100 

HALOE 

8 

80 

— 

NEPS 

— 

20 

— 


Table 13c: Subassembly Contributions to 

Cumulative Contaminant Deposition (Case #3) 


Subassembly 

Percent Coni 
GMW = 100 

.ribution to Conta 
GMW = 150 

minant Flux 
GMW = 200 

Solar Panel 
(ft = 0.412 steradians) 

75 

— 

94 

Instrument Module 
(ft = 0.0354 steradians) 

13 

— 

6 

HALOE 

12 

100 

— 


Table 14: Cumulative Volatile Organic Contaminant Deposition at HALOE Aperture 

(values based on do in parentheses) 


Species 

Cumulative Depth Over 35 Months [A] 


Case #1 

Case #2 

Case #3 

ficticious organic volatile species 

9496 

9787 

4663 

representing all GMW =100 materials 

(1602) 

(1652) 

(787) 

ficticious organic volatile species 

683 

408 

524 

representing GMW =150 material 

(173) 

(103) 

(133) 

ficticious organic volatile species 

13355 

489 

7460 

representing GMW = 200 material 

(4507) 

(165) 

(2518) 

total extrapolated 

23534 

10684 

12647 

accumulation per run 

(6282) 

(1920) 

(3438) 
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Table 15: Cumulative Volatile Organic Contaminant Deposition 
at HALOE Aperture — Arrival at Case- Averaged Total Deposition 

(values based on do in parentheses) 



Item 


Amount Based 

Amount Based 





On / = 10A [A] 

On do = 3.0A [A] 



Average for runs # 1 &; #2 

17109 

(4101) 



Avg. X (1 - Tgtow/T’orbit) 

2673 

(641) 



Average for run #3 

12647 

(3438) 



Avg. X Tgtow/^orbit 

10671 

(2901) 



Case-averaged Total Deposition 

13344 

(3542) 


1 

Table 16: Simulation Set-Up and Performance 


Item 

Inner Domain 

Outer Domain 

[sec] 

1.35x10-* 

2.69x10-* 

FNUM 

2.08xl0 13 

4.17xl0 13 

Volume [m 3 ] 

11.21x6.97x15.75 

21.20x21.20x21.20 

CCM Resolution 

74x46x 104 

70x70x70 

FCM Resolution w/in CCM 

3x3x3 

lxlxl 

No. of cells 

34,900 

7,530 

No. of simulated molecules 

329,000 

404,000 

Total CPU time [Ars] 

191 

150 

Memory required [Mbytes] 

60 

40 

Performance [psec/ particle/ At] 

44 

29 





92 NASA TM 109069 



Figure 1. UARS model used in DSMC simulation: (a) +Z side, Case 1; (b) -Z side, Case 1; (c) -Z side, 
Case 2; (d) -Z side, Case 3. Key: 1 -Cryogenic Limb Array Etalon Spectrometer (CLAES), 2-Solar panel 
(SP), 3- Multi -mission Modular Spacecraft (MMS), 4-Nadir Energetic Particle Spectrometer (NEPS), 5- 
Zenith Energetic Particle Spectrometer (ZEPS), 6-High Gain Antenna (HGA), 7-Solar /Stellar Positioning 
Platform (SSPP), and 8-Instrument Module (IM). 
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CPU #2 



Figure 2. Schematic representation of computational fiowfield volume used with parallelization routines. 



Figure 3. Representative cross-sections of computational grid for DSMC simulation; (a) inner domain at Z 
= 0, (b) inner domain at Y = 0. 



Figure 4. Representative cross-section of UARS geometric model discretized for incorporation into DSMC 
code. Outward-pointed normals included with surface pixel locations. 
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Figure 5. Typical behavior for number of particles per timestep versus simulated elapsed time as steady 
state is approached. 



F, S ure Typical distribution of particles per inner-domain cell before and after grid adaption. 
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Figure T. Case 1 number density contours of freest ream species 0, logarithmic scale nondimensionalized 
by freestream number density [0]^: (a) +Z side, (b) -Z side. 



I °9io( N n. / No,J 

lr 




Figure 8. Case 1 number density contours of neon, logarithmic scale nondimensionalized by [Oloc,: (a) -f Z 
side, (b) -Z side. 


K>9<N C o 5 /N 0 J 

R- 




MUN C02 /n 0 j 


Figure 9. Case 1 number density contours of carbon dioxide, logarithmic scale nondimensionalized bv 10} 
(a) +Z side, (b) -Z side. 
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Figure 10. Case 1 number density contours of GMW = 100 outgassing species, logarithmic scale nondimen- 
sionalized by [O]oo'- (a) +Z side, (b) - Z side. 
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Figure 11. Case 1 number density contours of GMW = 150 outgassing species, logarithmic scale nondimen- 
sionalized by [0]„: (a) +Z side, (b) -Z side. 
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Figure 12. Case 1 number density contours of GA4W ~ 200 outgassing species, logarithmic scale nondimen 
sionalized by [C)]*,: (a) +Z side, (b) ~Z side. 













Figure 13. Case 1 species number density contour maps emphasizing the HALOE aperture environment. 
Logarithmic scale nondimensionalized by freestream number density [O] (a) monatomic oxygen, (b) neon, 
(c) carbon dioxide, (d) GMW = 100 species, (e) GMW = 150 species, and (f) GMW — 200 species. 
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Figure 14. Case 2 species number density contour maps emphasizing the HALOE aperture environment 
Logarithmic scale nondimensionalized by freestream number density [Ojoo : (a) monatomic oxygen, (b) neon 
(c) carbon dioxide, (d) GMW = 100 species, (e) GMW = 150 species, and (f) GMW = 200 species. 
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Figure 15. Case 3 species number density contour maps emphasizing the HALOE aperture environment. 
Logarithmic scale nondimensionaJized by freestream number density [OJ^: (a) monatomic oxygen, (b) neon, 
(c) carbon dioxide, (d) GMW = 100 species, (e) GMW = 150 species, and (f) GMW = 200 species. 
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(b) 



(c) 


Figure 16. Number of surface collisions encountered by particles prior to striking the HALOE aperture: 
(a) Case 1, (b) Case 2, and (c) Case 3. 
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Figure 17. Case 1 histograms representing velocity distribution functions for particles striking the HALOE 
aperture: (a) X direction, (b) Y direction, (c) 2 direction. 
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Figure 18. Case 2 histograms representing velocity distribution functions for particles striking the HALOE 
aperture: (a) X direction, (b) Y direction, (c) Z direction. 
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Figure 19. Case 3 histograms representing velocity distribution functions for particles striking the HALOE 
aperture: (a) X direction, (b) Y direction, (c) Z direction. 
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(a) 



(b) 


Figure 20. Actual HALOE signal attenuation measurements between October 11, 1991 and October 20, 
1992 for (a) HF and (b) HCl species. 
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